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Abstract 

In this paper the problem of stabilizing large-scale systems by distributed 
controllers, where the controllers exchange information via a shared lim- 
ited communication medium is addressed. Event-triggered sampling schemes 
are proposed, where each system decides when to transmit new informa- 
tion across the network based on the crossing of some error thresholds. 
Stability of the interconnected large-scale system is inferred by applying 
a generalized small-gain theorem. Two variations of the event-triggered 
controllers which prevent the occurrence of the Zeno phenomenon are also 
discussed. 

1 Introduction 

We consider large-scale systems stabilized by distributed controllers, which com- 
municate over a limited shared medium. In this context it is of interest to reduce 
the communication load. An approach in this direction is event-triggered sam- 
pling, which attempts to send data only at "relevant times". In order to treat 
the large-scale case, input-to-state stability (ISS) small-gain results in the pres- 
ence of event-triggering decentralized controllers are presented. 
The stability (or stabilization) of large-scale interconnected systems is an impor- 
tant problem which has attracted much interest. In this context the small-gain 
theorem was extended to the interconnection of several £ p -stable subsystems. 
Early accounts of this approach are [28] (see also [22]) and references therein. 
For instance, in [28], Theorem 6.12, the influence of each subsystem on the 
others is measured via an £ p -gain, p£ [l,oo] and the £ p -stability of the inter- 
connected system holds provided that the spectral radius of the matrix of the 
gains is strictly less than unity. In other words, the stability of interconnected 
£ p -stable systems holds under a condition of weak coupling. 
In the nonlinear case a notion of robustness with respect to exogenous inputs 
is input-to-state stability (ISS) ([23]). If in a large-scale system each subsys- 
tem is ISS, then the influence between the subsystems is typically modeled via 
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nonlinear gain functions. Small-gain theorems have been developed for ISS sys- 
tems as well ([13, 14, 26]) and more recently they have been extended to the 
interconnection of several ISS subsystems ([7, 8]). For a recent comprehensive 
discussion about the literature on ISS small-gain results see [17]. 
In the literature on large-scale systems we have discussed so far, the communi- 
cation aspect does not play a role. If however, a shared communication medium 
leads to significant further restrictions, concepts like event-triggering become of 
interest. We speak of event-triggering if the occurrence of predefined events, as 
e.g. the violation of error bounds, triggers a communication attempt. Using this 
approach a decentralized way of stabilizing large-scale networked control systems 
which are finite £ p -stable has been proposed in [29, 31]. In these papers each 
subsystem broadcasts information when a locally-computed error signal exceeds 
a state-dependent threshold. Similar ideas are presented in [25, 30]. Numerical 
experiments e.g., [30] show that event-triggered stabilizing controllers can lead 
to less information transmission than standard sampled-data controllers. For 
consensus problems, event-triggered controllers are studied in [9]. 
One drawback of the proposed event-triggered sampling scheme is the need for 
constantly checking the validity of an inequality. A related approach which tries 
to overcome this issue is termed self-triggered sampling (see e.g., [3, 19]). 
From a more general perspective, the way in which the subsystems access the 
medium must be carefully designed. In this paper we do not discuss the problem 
of collision avoidance. This problem is addressed for instance in the literature 
on medium access protocols, such as the round-robin and the try-once-discard 
protocol. E.g., in [20] a large class of medium access protocols are treated 
as dynamical systems and the stability analysis in the presence of communi- 
cation constraints is carried out by including the protocols in the closed-loop 
system. This allows to give an estimate on the maximum allowable transfer 
interval (MATI) , that is the maximum interval of time between two consecutive 
transmissions which the system can tolerate without going into instability. The 
advantage of event-triggering lies in the possibility of reducing overall commu- 
nication load. However, if events occur simultaneously at several subsystems 
the problem of collision avoidance remains. We will discuss this in future work. 
The purpose of this paper is to explore event-triggered distributed controllers for 
systems which are given as an interconnection of a large number of ISS subsys- 
tems. Since input-to-state stability and finite L p stability are distinct properties 
for nonlinear systems, the class of systems under consideration in this paper dif- 
fers from the one in [29, 31]. Moreover, we use analytical tools which have been 
extended to deal with other classes of systems (such as integral-input-to-state 
stable systems [12] and hybrid systems [17]), and therefore the arguments in 
this paper are potentially applicable to a larger class of systems than the one 
actually considered here. 

We assume that the gains measuring the degree of interconnection satisfy a 
generalized small-gain condition. To simplify presentation, it is assumed fur- 
thermore that the graph modeling the interconnection structure is strongly con- 
nected. This assumption can be removed as in [8]. Since our event-triggered 
implementation of the control laws introduces disturbances into the system, the 
ISS small-gain results available in the literature are not applicable. An addi- 
tional condition is required for general nonlinear systems using event-triggering. 
This condition is explicitly given in the presented general small-gain theorem. 
Moreover, the functions which are needed to design the state-dependent trigger- 
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ing conditions are explicitly designed in such a way that the triggering events 
which supervise the broadcast by a subsystem only depend on local information. 
As an introductory example we explicitly discuss the special case of linear sys- 
tems, although for this class of systems the techniques of [29, 31] are applicable. 
As distributed event-triggered controllers can potentially require transmission 
times which accumulate in finite time, we also discuss two variations of the pro- 
posed small-gain event-triggered control laws which prevents the occurrence of 
the Zeno phenomenon. Related papers are also [10], [18]. 

Section 2 presents the class of system we focus our attention on, along with a 
number of preliminary notions and standing assumptions. The definition of the 
term event-triggered control can be found in Section 3. 

In Section 4 the notion of ISS-Lyapunov functions is presented. Based on this 
notion small-gain event-triggered distributed controllers are discussed in Sec- 
tion 5. The results are particularized to the case of linear systems in Section 2.1 
along with a few simulation results in Section 6. A nonlinear example together 
with simulation results is discussed in Section 7. 

The Zeno-free distributed event-triggered controllers are proposed in Section 8. 
The last section contains the conclusions of the paper. 

Notation No = N U {0}. R+ denotes the set of nonnegative real numbers, 
and R™ the nonnegative orthant, i.e. the set of all vectors of R™ which have all 
entries nonnegative. By || • || we denote the Euclidean norm of a vector or a 
matrix. 

A function a : R + — > R + is a class-/C function if it is continuous, strictly 
increasing and zero at zero. If it is additionally unbounded, i.e. lim r _5. +co a(r) = 
oo, then a is said to be a class-ZC^ function. We use the notation a € JC 
(a e /Ceo) to say that a is a class-/C (class-ZC^) function. The symbol JC U {0} 
(/Coo U{0}) refers to the set of functions which include all the class-/C (class-/Cco) 
functions and the function which is identically zero. A function a : M + — > M + 
is positive definite if a(r) = if and only if r = 0. We denote the right-hand 
limit by lkty x(t) — x(t + ). 

2 Preliminaries 

Consider the interconnection of N systems described by equations of the form: 



where i g JV := {1,2,..., N}, x = (xj . . . a;^) T , with x { e K" 1 , is the state 
vector and Ui € R m; is the zth control input. The vector e, with e = (ej . . . e^) T 
and € R ni , is an error affecting the state. We shall assume that the maps 
fi satisfy appropriate conditions which guarantee existence and uniqueness of 
solutions for inputs e. In particular, the /, are continuous. Also we assume 
that the gi are locally bounded, i.e. for each compact set K C R™ (n := Y^i=i n i) 
there exists a constant Ck with ||<7j(x)|| < Ck for each x € K. 
The interconnection of each system i with another system j is possible in two 
ways. One way is that the system j influences the dynamics of the system i 
directly, meaning that the state variable Xj appears non trivially in the function 
fi. The other way is that the controller i uses information from system j. In 



Xi 
Ui 



fi(x,Ui) 

gi(x + e) , 



(1) 
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this case, the state variable Xj appears non trivially in the function gi (and 
affects indirectly the dynamics of the system i). 

In this paper we adopt the notion of ISS-Lyapunov functions ([24]) to model 
the interconnection among the systems. 

Definition 1 A smooth function V : R" — > R + is called an ISS-Lyapunov 

function for system x = f(x, u) if there exist a\, ctz € /Coo and a.3, x e /C, such 
that for any x G R n 

ai(N|)<n*)<a 2 (N|) 
and the following implication holds for all x e K" and a// admissible u 

V(x) > X (\\u\\) => VV(x)f(x,u) < -a 3 (\\x\\) . 

It is well known that a system as in Definition 1 is ISS if and only if it admits an 
ISS-Lyapunov function. If there are more than one input present in the system, 
the question how to compare the influence of the different inputs arises. To 
answer this question we preliminary recall the notion of monotone aggregate 
functions from [8]: 1 

Definition 2 A continuous function fi : K™ — > R + is a monotone aggregation 
function if: 

(i) fj,(v) > for all and fj,(v) >0ifv>0; 

(ii) n{v) > ji(z) if v > z; 

(Hi) If \\v\\ — > 00 then ^l{v) — > 00. 

The space of monotone aggregate functions (MAFs in short) with domain R™ 
is denoted by MAF n . Moreover, it is said that [i e MAF™ if for each i = 
1,2,..., m, ^ e MAF n . 

Monotone aggregate function are used in the following assumption to specify the 
way in which systems are interconnected and how controllers use information 
about the other systems: 

Assumption 1 For i = 1,2, . . . , N , there exists a differentiable function Vi : 
W li — > ]R +7 and class-ICoo functions on\,a.ii such that 

MINI) < Vi( Xi ) < a^dl^iH) . 

Moreover there exist functions fa € MAF 2 at, Jij^ij € ^00 U {0} 7 di positive 
definite such that 

Vi(xi) >Mi(7a(^i(a;i)), ■ • • ,7iN(V N (x N )),Va(\\ei\\), ■ ■ ■ ,ViN{\\e N \\)) 
=>Wi(xi)fi(x,gi(x + e)) < -Oj(||a;j||) . 

1 In the definition below, for any pair of vectors v,z 6 I™, the notations v > z, v > z 
are used to express the property that Vi > Zi, Vi > Z{ for all i = 1,2,..., n. Moreover, the 
notation v ^ z indicates that v > w and v ^ w. 
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Loosely speaking, the function 7^ describes the overall influence of system j 
on the dynamics of system i, while the function 77^ describes the influence of 
the system j on the system i via the controller Qi. In particular, r]ij 7^ if and 
only if the controller Ui is using information from the system j. In this regard 
rjij describes the influence of the imperfect knowledge of the state of system j 
on system i caused by e.g., measurement noise. On the other hand, if i 7^ j 
and jij 7^ 0, then the system j influences the system i (either explicitly or 
implicitly). We assume that -fa = for any i. Observe that if the system i is 
not influenced by any other system j 7^ i, and there is no error on the state 
information Xi used in the control Uj, then the assumption amounts to saying 
that the system i is input-to-state stabilizable via state feedback. 

Remark 1 In general it is hard to design controllers that render the closed 
loop system ISS as we demand in Assumption 1. Though, there exist design 
techniques for special classes of systems. See e.g., [15, 16, 6] and the references 
therein. 

For future use we denote the set of states entering the dynamics of system i by 

= {j € -A/" : fi depends explicitly on Xj} , 

where explicit dependence of fi on Xj means that dfi/dxj ^ 0. Similarly for 
the controllers we denote 

= {j € M ■ g% depends explicitly on Xj} . 

It is also convenient to define the set of the controllers to which the state of 
system i is broadcast 

= {j e N '■ 9j depends explicitly on Xi} . 
2.1 The case of linear systems 

To get acquainted with the assumption above, we examine in the following 
example the case in which the systems are linear. 

Example 1 Consider the interconnection of N linear subsystems 

X} — ^^^=1 ^-ij^j ^illi 

Ui = ^2j = i Kij(xj + ej) . 

For each index i, we assume that the pairs (An,Bi) are stabilizable and we 
let the matrix Ku be such that An := An + BiKa is Hurwitz. Then for each 
Qi = Qj > there exists a matrix Pj = > such that AJ i P i + Pi An = —Qi 
leading to Lyapunov functions Vi(xi) — xj P t Xi. 
We consider now the expression Wi(xi)xi where 

Xi — y~^^-i (Ag + BiKij^jx j + BiKijCj 
='■ Sj=i AijXj + ^2j—i BijCj , 
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with Bij := BiKij and A^ := Aij + BiKij. 
Standard calculations lead to 



N 



N 



V^(x l )x J <-c l ||a; J || 2 +2||^||||P 4 ||( ^ ||^-||+^ ||B y | 

where 2 Cj = A m i n (Qj). Moreover, for any < q < the inequality 

iV AT 



E Il^-IIINI + EII^IIIN 



implies that 



The former inequality is implied by 



Vi( Xi )> \\Pi\ 



We conclude that (2) holds with 





2||P,|| 3/2 



£*i(r) = (ci-Cj)r 2 



[A„ m (P 3 )]!/2 

I Ik 



(3) 



It is important to remember that not all the functions 7^ and Tfo- are non-zero. 
Namely, 7^ (« 7^ j) is non-zero if and only if is a non-zero matrix. Similarly, 
77ij 7^ if and only if B^ 7^ 0. 



3 Event-triggered control 

In this paper we investigate event-triggered control schemes. Such schemes (or 
similar) have been studied in [3, 19, 25, 29, 30, 31]. 

We consider systems as defined in (1). Combined with a triggering scheme the 
setup under consideration has the form 

&i fi {<Ei Hi ) 

u, = g l (x + e) 

x = (4) 

e = x — x 

with triggering condition 

T 4 (x 2 , ei )>0. (5) 

2 For symmetric Qi we let A m i n (Qi) denote the smallest eigenvalue of Qi. 
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Here Xi is the state of system i <G Af, x is the information available at the 
controller and the controller error is e = x — x. We assume that the triggering 
function Tj are jointly continuous in a;,, ej and satisfy T t (xi, 0) < for all Xi ^ 0. 
Solutions to such a triggered feedback are defined as follows. We assume that 
the initial controller error is e = 0. Given an initial condition x we define 

ii := inf{i > : 3i e Af s.t. T^x^t), e^t)) > 0} . 

At time instant t\ the systems j for which Tj(xj,ej) = broadcast their respec- 
tive state Xj to all controllers with j e C(i). In particular, ej(tf) = for these 
indices j. 

Then inductively we set for k = 1, 2, . . . 

tfc+i := inf{* > t k : 3i € Af s.t. T,(x 4 (i), e 4 (t)) > 0} . 

We say that the triggering scheme induces Zeno behavior if for a given initial 
condition xq the event times tk converge to a finite t* . 

Remark 2 

• One of the proposed triggering schemes in this paper uses the information 
di which is an estimate of available at system i. For this scheme the 
triggering condition will be replaced by Ti(xi,ei,di) > 0. 

• The condition e = is used for simplicity. The triggering scheme uses 
implicitly that system i knows its state Xi and the error at the controller 

(and possibly the estimate di if this is used). It would therefore be 
sufficient to have an initial condition where system i is aware of e M = 
Xi(0) — Xi(0). However, such an assumption is most likely guaranteed by 
an initial broadcast of all states of the subsystems. But then e = is 
plausible. 

• It is a standing assumption in this paper that information transmission 
is reliable, so that broadcast information is received instantaneously and 
error free by the controllers. If this is not the case, additional techniques 
as studied e.g. in [27] have to be employed. This will be the topic of future 
research. 

• In many useful triggering conditions we have that 7^(0,0, df) — 0. If 
the system were to remain at x = this would lead to a continuum 
of triggering events, which do not provide information. To avoid this 
(academic) problem we propose to add the condition that information is 
broadcast once Xi reaches the state zero, but no further transmission by 
system i occurs as long as it stays at zero. 

• For simplicity, we assume x = in between triggering times. Usually, this 
is referred to as zero order hold. 

Other techniques are also possible, which could lower the triggering fre- 
quency. Consider for instance the case that each controller has a model 
for the dynamics of each other subsystem. Then each controller could use 
these models to calculate x rather than keeping it constant. Another ap- 
proach would be to extrapolate x linearly with the help of the last values 
for x. This is known as predictive first order hold. Both techniques would 
lead to x ^ 0. The considerations in this paper would also hold true for 
these cases with slight modifications of the proofs. 
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4 ISS Lyapunov functions for large-scale sys- 
tems 



In this section we review a general procedure for the construction of ISS Lya- 
punov functions. In particular, we extend recent results to a more general case 
that covers the case of event-triggered control. 

Condition (2) can be used to naturally build a graph which describes how 
the systems are interconnected. Let us introduce the matrix of functions V e 
(/Coo U{0}) JVxAr defined as 



r = 



/ 

721 



712 





713 
723 



\ 7jvi 7N2 7W3 



7ijv \ 

72 N 

J 



NxN 



Following [8], we associate to T the adjacency matrix A T = [a^] e {0,1} 
whose entry ay is zero if and only if 7^ = 0, otherwise it is equal to 1. Ay 
can be interpreted as the adjacency matrix of the graph which has a set N of 
TV nodes, each one of which is associated to a system of (1), and a set of edges 
£ C J\f x J\f with the property that (j, i) S £ if and only if a,j = 1. Recall that 
a graph is strongly connected if and only if the associated adjacency matrix is 
irreducible. In the present case, if the adjacency matrix Ap is irreducible, then 
we say that F is irreducible. In other words, the matrix of functions T is said 
to be irreducible if and only if the graph associated to it is strongly connected. 
For later use, given e MAFjv, Jij € /Coo U {0}, it is useful to introduce the 
map r M : 



1+ defined as 



r„(r) = 



/ /"i(7ii(n),--.,7uv(r-iv)) \ 



V MJv(7ivi(ri), . . . ,lNN(r N )) J 



Since the functions which describe the interconnection of the system are in 
general nonlinear, the topological property of graph connectivity may not be 
sufficient to ensure stability properties of the interconnected system. There 
must also be a way to quantify the degree of coupling of the systems. In this 
paper, this is done using the following notion: 

Definition 3 A map u <G JC^ is an Cl-path with respect to if: 

(i) for each i, the function a^ 1 is locally Lipschitz continuous on (0, 00); 

(ii) for every compact set K C (0, 00) there are constants < c < C such that 
for all i = 1,2, ... ,N and all points of differentiability of a' 1 we have: 

< c < {o-^Yir) < C , VreK; 



(Hi) T^(a(r)) < <r{r) for all r > 0. 

Condition (iii) in the definition above amounts to a small-gain condition for 
large-scale non-linear systems (in other words, condition (iii) requires the degree 
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of coupling among the different subsystems to be weak. For a more thorough 
discussion on condition (hi) see [8]). To familiarize with the condition, take the 
case N — 2 and \i\ = fi 2 = max (it is not difficult to see that the function 
maxi<i<jvri belongs to MAFn). Then 

r (r) = ( 7l2(r2) ^ 
^ (r >-{ 721 (n) ) ■ 

We want to show that there exists a <G K?^ such that T ll (a(s)) < <r(s) for all 
s > if and only if 712 o 721(7*) < r for all r > (the latter can be viewed as 
a small-gain condition for the interconnection of two ISS-subsystems). To this 
purpose, choose 




where 721 < cr 2 < 7 12 X - As a consequence of this choice, T fJ ,(a(s)) becomes: 

By construction, 712(0-2(5)) < s = 01 (s) and 721 (s) < er 2 (s), he. T ll (a(s)) < 
<j(s) for all s > 0. 

Strong connectivity of T and an additional condition implies a weak coupling 
among all the systems, in the following sense (see [8] for a proof and a more 
complete statement): 

Theorem 1 Let T e (/C^ U {0}) jVxjY and n e MAF^. IfT is irreducible and 
^ id 3 then there exists an Q-path a with respect to T^. 

Remark 3 In fact, the irreducibility condition on T is a purely technical as- 
sumption. A way how to relax it can be found in [8]. 

The small gain condition stated above is reformulated in the following assump- 
tion to take into account the case in which the error inputs are present in the 
system: 

Assumption 2 There exist an Q-path a with respect to T M and a map ip G 

(/Coo U {0} ) N x N such that: 

T ll (a(r)Mr))<v{r), Vr > , (6) 
where r jU (cr(r), <p(r)) is defined by 

(Mi (7ii(cti (r)), ■■,Jin(o- N (r)),fii(r), ..,(p 1N (r)) \ 
; 
Mjv(7ivi(o-i(r)), ..,jNN(o- N (r)),(pNi(r), ..,^jvjv(r)) / 

3 T M ^ id, means that for all s ^ r M (s) 2 s, i- e - f° r a H s S such that s ^0 there exists 
i 6 Af for which . . . , Sjv) < Si- 
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Remark 4 We remark that in the case ^ = max for each i E Af, one can 
exploit the degree of freedom given by ip in such a way that the condition (6) 
boils down to the small-gain condition r M (r) ^ r. In fact, once an fi-path has 
been determined, if the small-gain condition is true then it suffices to choose ifij 
such that, for any i, j e N, tpij < ^ik ° o k for some k e N. For a more general 
discussion on the fulfillment of (6) as a consequence of the small-gain condition, 
we refer the interested reader to [8], Corollaries 5.5-5.7. 

Remark 5 Observe that jij describes the influence of system j on the dynamics 
of system i either directly or through its controller gi. Hence for i ^ j the gains 
-fij ^ if and only if j € S(i) or j e C(i). Analogously, rjij ^ if and only 
if j e C(i), meaning that the controller i depends explicitely on the state of 
system j. Because the ipij from Assumption 2 describe the gains for the error 
input, there is no loss in generality if we set conventionally ip^ — if r]ij = 0. 
From the definition of C and Z it is evident that j € C(i) is equivalent to 

i e z(j). 

5 Main results 

In our first result it is shown that a Lyapunov function V and a set of decentral- 
ized conditions exist which guarantee that V decreases along the trajectories of 
the system: 

Theorem 2 Let Assumptions 1 and 2 hold. Let V{x) = ma,x ie ^f o-^ 1 (Vi(x i )) 
and, for each j e Af, define: 

Xj = o-j o % , with fjj = ^max (p^ 1 o VlJ . (7) 

Then there exist a positive definite a : E + — > E + such that the condition 
Vi(xi)>Xi(\\ei\\), yi E Af (8) 

implies 

(p,f(x,g(x + e))) < -a(\\x\\), Vp e dV{x) , 
where dV denotes the Clarke generalized gradient 4 and 

h{x,gi(x + e)) 



f(x,g(x + e)) 



f n {x,g n (x + e)) 



Proof: For each x, let Af(x) C Af be the set of indices i for which V(x) = 
o-^iV^Xi)). Let i e A/"(x) and set r = V(x). Then 

K(^) = Gi{r) > r^(o-(r),<p(r)) (Q) 
= l^i(jii(o-i{r)),..,jiN{o-N{r)),<Pii(r),..,ipi N {r)). 



4 We recall that by Rademacher's theorem the gradient VV of a locally Lipschitz function 
V exists almost everywhere. Let N be the set of measure zero where VV does not exist 
and let S be any measure zero subset of the state space where V lives. Then dV(x) = 
co{limi^ +00 SJV(xi) : Xi x, Xi N Xi S}. 
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Observe that by definition of V(x), for any i G N{x) and any j G J\f, 

liiWAr)) = H^j(V(x))) > ^(a^iVjix^)) = ^(V^ j) . (10) 

Note that for j ^ C(i) we have <pij — and r)ij = 0. Hence for j £ C(i) it holds 
trivially that 

^•M > %(IMI). (ii) 

This is also true if j G C(i) (or cquivalcntly i G Z(j)). In fact, since for any 

Vj(xj) > Xj{\\ej\\) , Xj=Vj°Vj 
we have, using the definition of V, (8) and (7), that 

<Pij(r) = <Pij{V(x)) > ifijia^^jixj))) > fijicrj 1 o CT j (%(||e J -||))) 

> ^(crT 1 off^ipr. 1 ojjydlejH))) = »fa(ll e jll) • (I 2 ) 

Observe that ^i(w) > Pi(.z) for all w > z G since ^ G MAF 2 n an d as a 
consequence of Definition 2, (ii). Since r = V(x) > a~ 1 (Vi(x i )) for all i G N ', 
by (10), (11) and (12), 

Mi(7ii(o"i(r)), . . . ^ iN {^N (r)), ViiWi • • • i¥>»Jv( r )) > 

/^(TuOMzi)), • • - ,7»Ar(Viv(a;jv)),»?ii(||ei||),.. . ,»?ijv(||ejv||)) • (13) 

The inequality above and (9) yield that for each i G Af(x) 

Vi(xi) > /ii(7ii(ai(r)), .., 7iJ\r(o-jv(r)), p»i(r), ..,<p iN (r)) 

> ^(iti(Vi(xi)), . . .,7ijv(Viv(a;jv)),»?ii(||ei||),.. . ,%jv(||ejv||)) • (14) 

Hence, by (2), 

VK(ii)/j(a,ft(!C + e)) < -a t (\\x t \\) 

for all i G Af(a;). 

We now provide a bound to (p, fi{x,gi(x + e))) for each p G dcr," 1 ^^)) and 
i G Af(x). Observe that er _1 is only locally Lipschitz and the Clarke generalized 
gradient must be used for a^ 1 (V i (x i )). Fix Xi and let p > be such that 
= p. Define the compact set K p = {Vi(xi) G M+ : p/2 < \\xi\\ < 2p}, and 

let 

c p = mm (err 1 )' W . C P = mfK" 1 )' (r) , 

where c p > by definition of the J7-path er. Bearing in mind that ||a;i|| = p, for 
each p G 9<7,~ (V^(o;j)) there exists 7 P G [c p ,C p ] such that p = j p Wi(xi), and 
(p,fi{x,gi(x + e))) = 7 P V^(xi) • fi(x,gi(x + e)) < ~7 P a 4 (p) < -c p a 4 (p). 
Set 5,(p) := c p ai(p), which is a positive function for all positive p. Also set 

a(r) := min{a;(||a;j||) : r = ||x|| , i G A/"(a;)} . 

Then, for each p G do-r^V^)), (p, fi(x,gi(x + e))) < — cti(j |a? 4 ||) < -a(||a;||). 
This in turn implies ([8]) that for eachpG <9F(x) (p, f(x, g{x + e))) < —a{\\x\\). 

□ 
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In the rest of the section we discuss an event-triggered control scheme for the 
system (4) with triggering conditions that ensure that the condition on the state 
x and the error e as in Theorem 2 are satisfied. 



Theorem 3 Let Assumptions 1 and 2 hold. Consider the interconnected system 

Xi(t) = fi{x(t), gi (x{t))) , i e Af , (15) 

as in (4) with triggering conditions given by 

Ti{xi,ei) = Xi(INI) - Vi{xi) , 

with Xi defined in (7) for all i G Af. Assume that no Zeno behavior is induced, 
i. e. the sequence of times tk , where the tk 's are defined by the triggering con- 
ditions Tj as discussed in Section 3, has no accumulation point or is a finite 
sequence for all i G Af. Then the origin is a globally uniformly asymptotically 
stable equilibrium for (15). 

Proof: To analyze the event-based control scheme introduced above, we 
define the time-varying map f(t, x) = f{x, g(x + e(t))). The map f(t, x) satisfies 
the Caratheodory conditions for the existence of solutions (see e.g. [4], Section 
1.1). Because of the conditions on / (sec Section 2), the solution exists and 
is unique. Along the solutions of x = f(t,x), the locally Lipschitz positive 
definite and radially unbounded Lyapunov function V(x) introduced in Theorem 
2 satisfies 

t" 



V{x{t"))-V{x{t'))= jV{x{t))dt 



for each pair of times t" > t' belonging to the interval of existence of the solution. 
Moreover, by a property of the Clarke generalized gradient ([5], Section 2.3, 
Proposition 4), for almost all t G R+, there exists p G dV(x{t)) such that: 

jV(x(t)) = {p,f(t,x(t))). 

Note that the triggering conditions Ti{xi,e{) — X«(ll e i|l) — Vi( x i) > ensures 
that Vi(xi) > Xi(ll e i|l) f° r an positive times. Hence we can use Theorem 2 
together with the definition of f(t, x), to infer (see [21], Section IV. B, for similar 
arguments) 

t" 

V{x{t"))-V{x{t'))<- a{\\x{t)\\)dt . 

We can now apply [4], Theorem 3.2, to conclude that the origin of x = f(t,x), 
and therefore of x = f{x,g{x)), is uniformly globally asymptotically stable. □ 
The assumption that no Zeno behavior is induced is quite strong. One possibility 
is to cast the event-triggering approach in the framework of hybrid systems and 
study the asymptotic stability of the system in the presence of Zeno behavior (see 
[11], pp. 72-73 for a discussion in that respect). Another possibility is to extend 
the solution ([2]). Let t* be the accumulation time such that linifc^ +00 tk = t* . 
Since the Lyapunov function is decreasing along the solution over the interval 
of time [0, t*), then limfc^ +00 V(x(tk)) exists and is finite. Let us denote this 
limit value as V* . If V* ^ 0, then one can pick a state x* such that V(x*) = V* 
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and consider the solution to the system (15) with initial condition x* . If Zeno 
behavior appears again, one can repeat indefinitely the same argument and 
conclude that V(x(t)) converges to zero either in finite or in infinite time, with 
x{t) obtained by the repeated extension of the solution after the Zeno times. 
However, this approach may raise a few implementation issues, such as the 
detection of the Zeno time and the choice of the new initial condition x* at the 
Zeno time, and may discourage to follow this path. For this reason, slightly 
different triggering conditions which rule out the possibility of Zeno behavior 
are introduced in Section 8. 



6 An example 

Consider the interconnection of linear systems as in Section 2.1 

JV JV 

±i = ^ AijXj + X! Bijej i£jV, (16) 

with An Hurwitz for i G J\f. In order to apply our event-triggered sampling 
scheme, we first have to check the conditions of Theorem 2. As verified in 
Section 2.1, Assumption 1 holds for system (16) with each Lyapunov function 
given by Vi(xi) = xj P i x l . 

To check Assumption 2 we recall Lemma 7.2 from [8]: 

Lemma 1 Let a e fCoo satisfy a(ab) = a(a)a(b) for all a, b > 0. Let D = 
diag(a), G e R nxn , and be given by 

T,(s) = D- 1 (GD(s)). 

Then ^ id if and only if the spectral radius of G is less than one. 

It is easy to see that for the linear case from Section 4 with entries from (3) 
is of the form of Lemma 1 with a(r) := y/r and 

G .. = g Il^-H nine* 

13 ~ ~ Cl [A min (P,)]V 2 ' ^3,hJ€M 

and zeros as diagonal entries. In other words, Jij(r) = Gija(r). Let us assume 
that the spectral radius of G is less than one. For the case of linear systems an 
f2-path is given by a half line in the direction of an eigenvector s* of a matrix G* 
which is a perturbed version of G (for details, see the proof of [8, Lemma 7.12]). 
Denote this half line by a(r) := s*r. 
To show a way to construct a ip for which 

T^a(r),<p(r)) < a(r) , Vr > (17) 

holds, consider the ith row of (17) and exploit the fact that the fi-path is linear: 

A 2||P i || 3 / 2 p„|| /— " . \ t 
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Bearing in mind that ^ if and only if j £ E(i) or j £ C(i) (see 

last paragraph of Example 2.1 together with the definition of the set X), the 
inequality can be rewritten as 



E 

vje(s(i)uc(i))\< 



2||fill 3/2 ll^ll 

Ci [XminiPj)} 1 / 2 



N 



H*+E^m < 



(18) 



If we make the choice <fij{r) — a%j-\/f for all j € C{i) and <fij(r) = 
otherwise, we obtain 



E 

j'ec(*) 



E 



2||^ll 3/2 HA 



je(2(i)uc(t))V 



5* =: ft. (19) 



It is worth noting that pi > by the spectral condition on G. 
Assume without loss of generality that (E(z) U C(i)) \i ^ (if not, (18) trivially 
holds). Note that it would be sufficient to assume irreducibility of the intercon- 
nection structure to ensure U C(i)) \ i ^ 0. 

Without further knowledge of the system, we choose for j £ C(i) for the gains 
tpij(r) := \ where \C(i)\ denotes the cardinality of the set C(i), to ensure 
that (17) holds. If j C(i) set tpy = 0. Simulations suggest that it might be 
better to not choose the uniformly, but to relate them to the system matrices 
(in particular, to the spectral radii of the coupling matrices Bij — BiKij). 
Now we can calculate the trigger functions Xi as m Theorem 2 by using the 
f2-path and the (fij from above. The map fji is calculated using the 77^ from 
(3). Stability of the interconnected system is then inferred by Theorem 3. 
To illustrate the feasibility of our approach we simulated the interconnection 





Figure 1: Trajectories of the intercon- Figure 2: Trajectories of the intercon- 
nected system with periodic sampling nected system with event-triggering 

of three linear systems of dimension three. The entries of the system matrices 
are drawn randomly from a uniform distribution on the open interval (—5,5). 
We repeat this procedure until the spectral radius of the corresponding matrix 
G is less than one. 

In Figure 1 new information is sampled every three units of time. Which sys- 
tem has to transmit information is decided by a round robin protocol (i.e., first 
system one, than system two, system three and again system one and so on). 
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Figure 3: 33 periodic (red dots) and 22 (blue stars) events at the beginning of 
the simulation 

3*- • * • x- • • * • x* • * • + • * • 

u 

"£2" + • -X • •* • *^ • • X- • X • • X- • 

M 

■| I— « !*• X— ± » : — X • : • 1 X* 1 • 1 — • 1 

1500 1510 1520 1530 1540 1550 1560 1570 1580 1590 t 1600 



Figure 4: 34 (red dots) periodic and 19 (blue stars) events in the middle of the 
simulation 



X- • • X- 



> -X • X- 



• X- • x» 



3 

i 2 

1 

2800 2810 2820 2830 2840 2850 2860 2870 2880 2890 t 290° 



Figure 5: 34 periodic (red dots) and 21 (blue stars) events at the end of the 
simulation 



In Figure 2 our event-triggered sampling scheme is used. 

Over the range of 3000 units of time the system with periodic sampling trans- 
mitted 1000 new information, whereas in our scheme the events were triggered 
only 595 times. By looking at Figure 1 and Figure 2 it seems like the periodic 
sampled system converges a bit faster. Indeed, the systems state norm of the 
event-triggered system at time 3000 is already reached by the periodic sampled 
system after 2486 (i.e., 828 periodic samplings) units of time. But still the 
number of triggered events (595) is smaller than the number of periodic events 
(828). 

A representation of how the the different systems (1,2, or 3) sample their state 
is depicted in Figures 3-5. The first picture shows the sampling behavior at the 
beginning (t 6 [0, 100]) of the simulation. The other two are from the middle 
(t e [1500, 1600]) and the end (t E [2800,2900]) of the simulation, respectively. 
There is a small overshoot for some of the trajectories in Figure 1. This behav- 
ior cannot be seen in Figure 2, because in the event-triggered implementation 
information is transmitted more frequently at the beginning by systems 2 and 3, 
whereas in the periodic implementation transmission starts (for systems 2 and 
3) a few samples later, as can be seen in Figure 3. The reason for this is that we 
set x = instead of initializing eo = 0. This is another possibility of initializing 
the controller (and hence the initial error) than the one described in Remark 2. 



7 A nonlinear example 

The following interconnection of N — 2 subsystems 



±1 = X1X2 + x\u\ 
±2 = x\+ u 2 
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is considered under the assumption that each controller can only access the state 
of the system it controls. The control laws are chosen accordingly as 

ui = -(xi + ei) , u 2 = -k(x 2 + e 2 ) , k > . 

Let Vi(xi) = \x 2 for i = 1,2. Then 

Vxixx) := Wi(xi)(-a;? + x l2 ; 2 - xfei) < a;? (--a;? + |z 2 | + ^i) ( 20 ) 

from which we can deduce 

^?>N + ^e? => V^i(a;i)<-^. 

Since the left-hand side of the implication is in turn implied by V\{x\) > 
max{-y/32V2(x2), 2ef }, this shows that the first subsystem fulfills Assumption 1 
with 



= max , 7n (r) = , 712(7") = VMr , r?n(r) = 2r 2 , r?i 2 (r) = . (21) 
Similarly 

V 2 (x 2 ) := W 2 {x 2 ){x\ - kx 2 - ke 2 ) < |a; 2 |(-fc|a;2| + x\ + k\e 2 \) 
and therefore 

32 • k 

V 2 (x 2 ) > max{plf(a:i),8ei} => V 2 (.t 2 ) < --a| , 

i.e. the second subsystem satisfies Assumption 1 with 
32 

/i 2 = max , 721 (r) = ^r 2 , 7 22 (r) = , ?? 2 i(r) = , ry 22 (r) = 8r 2 . (22) 

As discussed in Section 4, in the case of N = 2 the f2-path can be chosen as 
a 1 = Id and 7 21 < a 2 < 7 r2 1 . Provided that fc > 32, one can set a 2 (r) = a 2 r 2 , 
with a 2 e (p, jg). If y € (/Coo U {0}) JVxAr is additionally chosen as 

32 

ipu(r) = V32ar , tp 12 = <p 21 = , <^ 22 (r) = ^r 2 , 

then Assumption 2 is satisfied. In view of the choice of a, [i and ip, the require- 
ment (6) boils down to the condition T^(cr(r)) < u{r) which is equivalent to the 
small-gain condition 712 o 7 21 < Id (see Section 4). This small-gain condition 
is fulfilled by the choice of k, since 712 o 7 2 i(r) = j^r. Hence Theorem 2 ap- 
plies and provides an expression for the functions \i used in the event-triggered 
implementation of the control laws. The functions are given explicitly by 

1 7f 2 h 2 
Xi (0 := ~^=_r , X2(r) := —^-r . 

Simulation results for the initial condition 27(0) = —4, x 2 (0) = 3, 27 (0) = —4 
and i 2 (0) = 3 can be found for t € [0, 0.5] in Figures 6-8. 

The trajectory of the first system is given in blue and for the second system in 
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green. The input is calculated using the red and turquoise values accordingly. 
Figure 8 shows the event triggering scheme from Theorem 3. After 0.5 seconds 
12 events were triggered. Recall that we start with e(0) = 0. This initialization 
is not counted as an event. The shortest time between two events is 0.0155 
seconds. In Figure 6 we used a periodic sampling scheme with a sampling 
period of 0.0155 seconds leading to a total of 66 samples. Compared to the 
event triggering scheme no major improvement in the performance of the closed 
loop system can be recognized. Although in the periodically sampled system 
more than five times the amount of information was transmitted. 
If we use a Round Robin protocol with 12 periodic samples as in Figure 7 the 
sampling scheme is not able to stabilize the system. In Figure 9 the Lyapunov 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 . 0.5 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 , 0.5 



Figure 6: Trajectories of the intercon- Figure 7: Trajectories of the intercon- 
nected system with periodic sampling nected system with periodic sampling 
(66 samples) (12 samples Round Robin) 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 . 0.5 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 , 0.5 



Figure 8: Trajectories of the inter- Figure 9: Lyapunov function and 
connected system with event-triggered X2 (|| ^2 1| ) for the second subsystem 
sampling (12 events) 

function for the second subsystem together with X2(|| e 2||) is plotted. Every time 
the red curve (the error function) hits the blue line, the error is reset to zero. 

8 On Zeno-free distributed event-triggered con- 
trol 

The aim of this section is to show that it is possible to design distributed event- 
triggered control schemes for which the accumulation of the sampling times in 
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finite time does not occur. The focus is again on the system (1), namely: 



Xi = fi(x,gi(x + e)) . 



(23) 



We present two different approaches for a Zeno free event-triggered control 
scheme. The first is based on a practical ISS-Lyapunov assumption whereas 
the second tries to lower the amount of event-triggering by reducing unneeded 
information transmissions. 

8.1 Practical Stabilization 

Here we adopt a slight variation of the input-to-state stability property from 
Assumption 1. 

Assumption 3 For i € Af, there exist a differentiable function Vi : R ni — > R + , 
and class-K-oo functions an^a^ such that 



Moreover there exist functions /x, € MAF2JV; <= ^oo, for j € Af, positive 

definite functions a, and positive constants a, for i G Af, such that 



Vi(xi) > max{^(7 ll (V r 1 (x 1 )),. . . , 7 W (Vjv(ijv)), %i(||ei||), ■ • • , %jv(| |ejv| I)), cj 



Remark 6 Systems satisfying Assumption 3 arc usually referred to as input- 
to-state practically stable (ISpS) ([14]). 

We now state a new version of Theorem 2 for system (23) . 

Theorem 4 Let Assumptions 2 and 3 hold. Let V(x) = max ie ^f a^ 1 (Vi(x i )). 
Assume that for each j £ Af , 



MINI) < Vi(xi) <Oi2(IN|) . 



=> VV i (x i )f i (x,g l {x + e)) < -a^W) . (24) 



max{ a j 1 (V j (x j )),c j } > fjj{\\ej\\) , 



(25) 



where 




(26) 



Then there exists a positive definite a : R + — > R + such that 



(p, f(x, g(x + e))) < -a(\\x\\), Vp e dV(x) 



for all x = (xj xj . . . 



xJ[) T e {x : V(x) > c := maxjjc,,^ (cj)}}, where 
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Proof: Let Af(x) C A/" be the indices i such that V(x) = a i 1 (Vi(xi)). 
Take any pair of indices i,j € A/". By definition, V(x) > o-J 1 (Vj(xj)) and 

7«(^(^)))>7y(^(^))- (27) 
Let i e N{x). Then by Assumption 2, we have: 

Vi( Xi ) = <n(V(x)) > 

w(7iiM^(aO)), ■ ■ • , 7iJv(<M^0»0))> M^O*)), . . . , ^v(^))) • (28) 

Bearing in mind (27), we also have 

V i {x i )=u i {V{x))> 

Hi{li\{V\{xi)), . . . ,-f iN (V N {x N ))^ ll {V(x)) 1 . . . ,ip iN (V(x))) . (29) 

Let us partition the set Af := V U Q. The set V consists of all the indices 
i for which the first part of the maximum in condition (25) holds, i.e. i G 
V :<^> a^iViixi)) > q; also Q := N \ V. For all j e V we have by (25) 
a j~ (Vj( x j)) — Vj(\\ e j\\) an d hence using (26) (as mentioned before, the case 
j £ C(i) is trivial) 

<Pij(V(x)) > <pij o crJ^Vjixj)) > Vi j o %(|| ej ||) > 

<Pij ° Vy 1 ° %(IMI) = »?«(l|ej||) • (30) 
Assume now that V(x) > c. For all j e Q we have by (25) > %(||ej||) and so 
<Pij(V{x)) > <pij{c) > <pij{cj) > <Pij o %(|| ej ||) > TfadleJ). (31) 
Combining (30) and (31) we get for all j £ Af 

<Pij{V(x)) > f?ij(||ej||) , 

provided that (25) holds and that V(x) > c. Substituting the latter in (29) 
yields 

V i {x i )=a i {V(x)) > 

Vii-faiViixx)), . . .,7»Ar(Viv(a;jv)),»7ii(||ei||),.. ■ ,ViN(\\e N \\)) . (32) 

Because i € Af(x) and V(x) >c~ max^jci, a^ 1 (ci)}, we have 
V(x) = a^iViixi)) >c> o-- 1 ^) and finally we conclude 

Vi(xi) > a . 

The latter together with (32) is the left-hand side of the implication (24). Hence, 
Wi(xi)fi(x, gi(x + e)) < — q;^ (| jar^ | J ) for all i € Af(x). Now we can repeat the 
same arguments of the last part of the proof of Theorem 2, and conclude that for 
all x such that V(x) > c and for all p G dV(x), (p,f(x,g(x + e))) < — a(||a;||). 

□ 

Now that we established an analog to Theorem 2 for the case of practical sta- 
bility, we are able to infer stability of the closed loop event-triggered control 
system (4) without the assumption that Zeno behavior does not occur. 
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Theorem 5 Let Assumptions 2 and 3 hold. Consider the interconnected system 

x i (t) = f i (x(t),g i (x(t))), zeAA, (33) 

as in (4) with triggering conditions given by 

Ti{x t ,ei) = fn(\\ei\\) - max-ja^ 1 o Vi(xi),Ci} , (34) 

with fji defined in (26) for all i G N . Then the origin is a globally uniformly 
practically stable equilibrium for (33). 

Proof: Here we want to adopt the same line of reasoning as in the proof of 
Theorem 3. To this end, we have to make sure that by the triggering conditions 
given by (34) no Zeno behavior is induced. Note that in between triggering 
events e(t) = —x(t) for all t <E (tk,tk+i) by the definition of (4). The trig- 
gering conditions Ti(xi,ei) = fj%(\\ei\\) — maxjer" 1 o Vi(xi),di} > ensure that 
maxjer^ 1 o Vi(xi),&i} > %(||ei||) for all positive times. Following the same 
reasoning of Theorem 3, with the exception that we have to replace the applica- 
tion of Theorem 2 with the application of Theorem 4, one proves that V(x{t)) 
is decreasing along the solution x(t) on its domain of definition. Hence, x(t) 
is bounded on its domain of definition. Since maxjer^ 1 o Vi(xi(t)),di(t)} > 
j?i(||ej(t)||), then also e(t) is bounded and so is x (t) — x(t) + e(t). As ej(t^) = 
for each index j which triggered an event and e(t) is bounded in between events 
(e(t) = —x(t) = —f(x(t),g(x(t)))), the time when the next event will be trig- 
gered by system j is bounded away from zero because the time it takes ej to 
evolve from zero to Cj is bounded away from zero. Hence, either there is a finite 
number of times tk or tk — > oo as k goes to infinity. The solution x(t) is then 
defined for all positive times and Theorems 3 and 4 allow us to conclude. □ 
In hybrid systems, the practice of avoiding Zeno effects while retaining stability 
in the practical sense is referred to as temporal regularization (see [11], p. 73, and 
references therein). Here, the regularization is achieved via a notion of practical 
ISS. In the context of event-triggered £2-disturbance attenuation control for 
linear systems temporal regularization is studied in [10]. 



8.2 Parsimonious Triggering 

In Section 8.1 a way to exclude the occurrence of Zeno behavior for the price of 
practical stability rather than asymptotic stability was shown. Here we want to 
provide a way to exclude the Zeno effect by introducing a new triggering scheme, 
but still achieving asymptotic stability. The main idea behind the new triggering 
scheme, which will be introduced in Theorem 7 is that if the error of the zth 
subsystem is bigger than its Lyapunov function but still small compared to the 
Lyapunov function of the overall system, no transmission of the ith subsystem 
is needed. 

For future use we need also a slight variation of Theorem 2. Here we exploit 
the fact that we can either compare each state to its corresponding error (as 
in Theorem 2) or each error to the Lyapunov function of the overall system as 
shown in the next theorem. 

Theorem 6 Let Assumptions 1 and 2 hold. Let V(x) = ma,x ie ^f cr^ 1 (Vi(xi)) 
and, for each j e J\f, define: 

Vj = max ip^ 1 o Vij . (35) 
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Then there exist a positive definite a : R + — » R+ such that the condition 

V{x) > fjj(\\ej\\), Vj e TV (36) 

implies 

(p,f(x,g(x + e))) < -a(\\x\\), Vp e dV(x) . 

Proof: The proof follows by a slight modification of the proof of Theorem 2. 
For each x, let Af(x) C M be set of indices for which V"(x) = a," (Vi(xj)). 
It is sufficient to show that for all i G M{x), j e /V we have <fiij(r) > Vij(\\ e j\\)> 
with r = V~(x). 

First recall that for j C(i) the latter inequality trivially holds. So assume 
that j e C(i). Using (35) and (36), we have 

<Pij(V(x)) > <Pij(T)j(\\ej\\)) > <Pij o if- 1 o rjij ( 1 1 6j 1 1 ) = »*tf(||ej||) . 

i.e. (12) in the proof of Theorem 2. Then the argument after (12) can be 
repeated word by word. By previous arguments this concludes the proof. □ 
A triggering condition for the jth subsystem which yield the validity of condi- 
tion (36) would make the knowledge of the Lyapunov function V of the overall 
system to system j necessary. This would contradict our wish for a decentral- 
ized approach. 

The next lemma provides a decentralized way to ensure that condition (36) 
holds. To this end, we give an approximation of the other states (the W) which 
will be compared to the error instead of the Lyapunov function of the overall 
system. Appropriately scaled, IF is a lower bound on the Lyapunov function 
of the overall system and hence can be used to check the validity of (36). The 
important advantage is, that this approximation can be calculated by using only 
local information. Before we state the next lemma, define 

£?' x ' '■— (£i~> • • • >£j-ii x J ■ ■ ■ >£iv) 

as the vector £ where the jth component is replaced by Xj. The proofs of 
Lemma 2, 3, and 4 are postponed to the Appendix. 

Lemma 2 Let Assumptions 1 and 2 hold. Let V(x) = max ieA f a^ 1 (V i (x i )). Let 
dj e R + be an approximation of \\fj(x, gj(x + e))\\ with \ \\fj(x, gj(x+e))\\ — dj\ < 

Kj\\Xj\\. 

Assume that for j £ J\f there exist functions Qj : R N — > R such that V(x) > 
Vi(\\ e i\\) f or a M i 7^ j implies 

Q 3 (\\x l l... 1 \\x N \\)>\\f 3 (x,g J (x + e))\\. (37) 

Define 

W{j,Xj,dj) = min{max||£;|| : £ e A(j, Xj, dj)} 

with 

A(j, Xj ,dj) = {e*i : e j (U 1 \\,...,U N \\)> dj - ~ Kj \\xj\\} . (38) 
Then the conditions 

W(j, Xj ,dj) > ^J 1 o ^.(IMI) and aj\Vj{xj)) < Vj(\\ej\\), (39) 

with tpj = maxj^j ct" 1 o an imply 

V(x)>rjj(\\ej\\). 
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Remark 7 Lemma 2 presents a way of approximating the norm of the other 
states influencing the dynamics of a single subsystem. To this end an approx- 
imation of the derivative is used. Another possibility for achieving this goal 
would be to construct an observer, which gives an approximation of the inputs 
(the other states) by observing the dynamics of a single subsystem. 

Before we can state another event-triggering scheme, which does not induce Zeno 
behavior we have to formulate the observation that if Zeno behavior occurs, one 
of the states has to approach the equilibrium. 

Lemma 3 Consider a large scale system with triggered control of the form (4) 
satisfying Assumptions 1 and 2. Let \i, i € M be given by (7). Consider the 
triggering conditions given by 

Ti(xi,ei) = Xi(IMI) - Vi(xi) . 

If for a given initial condition x the triggering scheme Tj. induces Zeno behav- 
ior, then the corresponding solution 

Xi*(tk) -> • 

The next lemma provides an inequality for the state and the corresponding dy- 
namics. Besides the rather technical nature of Lemma 4, together with Lemma 3 
it forms the basis to be able to compare the ith state to the rest of the states 
as will be seen in Theorem 7. 

Lemma 4 Consider system 

x = f(x,g(x + e)) (40) 

as in (4) . If there are triggering instances tk — > t* for k — > oo and an index i 
such that Xi(tk) — > 0, then for all M > there exists a k* e N such that for 
some k > k* 

ll^fa +1 )-^fa)ll >M|kifa+i)|| , 
tfc+i — tfc 

It is of interest to note the following immediate corollary. 

Corollary 1 Under the conditions of Lemmas 3 and 4, assume that the func- 
tions @j from Lemma 2 satisfying (37) may be chosen to be Lipschitz and so 
that Oj(0, ...,0) = holds. Consider an initial condition x(0) = xq =/= 0. If 
there is Zeno behavior at t* , i.e. if there are triggering instances tk — > t* , then 
for the overall state x of (40) 

x(tfc) 7^ as k -> oo . (41) 

Proof: We first exclude that there is a s* <G [0, t*) such that x(s*) = 
0. Otherwise choose Lipschitz constants Li for Oj valid on the compact set 
{x(s) ; s € [0, s*]} and note that we have for each i almost everywhere on [0, s*] 

ll*<(*)ll = Il/i(^(*),ff*(a5(*) + e(*)))|| 

<e i (\\x 1 (t)\\,...,\\x N (t)\\)<L i \\ x ( t )\\ . (42) 
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Note that we can use Gi as a bound for the dynamics as in (37), because the 

validity of Vi{x{) > Xi(ll e i||) f° r au i trivially implies V{x) > j?i(||ej||). 

As (42) is true for all i, this implies ||x(i)|| < L||x(£)|| for L sufficiently large 

and almost all t e [0,s*]. It follows that ||x(s*)|| > e- Ls " \\x(0)\\ > 0, so that 

x(s*)^0. 

If x(tk) — > 0, then x(t) — > for t /• t* . Hence for each i, and k sufficiently large 
we have that (42) holds almost everywhere on (tk ,t*). As in the first part of 
the proof it follows that \\x(t*)\\ > e~ L (**- tfc )||:r(i fe )|| > 0, because by the first 
step of the proof x(tk) ^ 0. This contradicts the assumption that x(tk) — > 0. □ 
The rest of this section is devoted to constructing an event-triggered control 
scheme which ensures that Condition (36) holds. 

From Lemma 3 we know that if Zeno behavior occurs, then one of the subsys- 
tems approaches the origin in finite time. Corollary 1 shows that under certain 
regularity assumptions, a number of subsystems do not converge to as we ap- 
proach the Zeno point. Hence, from a certain time on, the Lyapunov function 
corresponding to the subsystem which tends to the origin does not contribute to 
the Lyapunov function for the overall system. As a consequence no information 
transfer from this subsystem is necessary using parsimonious triggering. This 
observation is made rigorous in the rest of the section. 

In the next theorem we use the triggering condition as in Theorem 3 but we 
add another triggering condition T i2 , which checks whether the ith subsystem 
contributes to the Lyapunov function of the overall system. It does so by com- 
paring the local error of system i with the approximation W of the other states 
as described in Lemma 2. The main idea is that if the dynamics of the ith 
system is large compared to its own state, other states must be large. As the 
correct value of the dynamics is not known to system i, an approximation of 
\\ii\\ is used. 

As the aim is to use only local information, we will use the difference quotients 
to approximate the size of the derivative at the triggering points. Furthermore, 
we do not wish to assume that all subsystems are aware of all triggering events. 
Hence in the following we will use the notation t\ to denote those triggering 
events initiated by system i. We define 

m = \mz^A=M (43) 

* *fe-l 

as the difference quotient approximating ||ii(t)|| after the triggering event t l kl . 
Adding the new triggering condition that uses (43) allows us to exclude the 
occurrence of Zeno behavior. A discussion about the new triggering condition 
can be found in Remark 8 and 9. 

Theorem 7 Consider a large scale system with triggered control of the form 
(4) satisfying Assumptions 1 and 2. Let V{x) — max^g^ a^ 1 (Vi(x i )). Define 

T il (x l ,e i ) = Xi(INI) - Vi(xi) 

with Xi as in Theorem 3 and 

T l2 {x t ,e l ,d l ) = ipr 1 o % ( 1 1 1 1 ) - W{i,x l ,d l ) , 

where tpi, W(i, x%, di) and fji are defined as in Lemma 2. Furthermore, assume 
that for alii G M the Gj from Lemma 2 and ofji are Lipschitz with Lipschitz 
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constant Li respectively Ki and that 0,(0, . . . , 0) = holds. 
Consider the interconnected system 

x i (t) = f i {x{t),g i (x{t))), i eJ\f , (44) 

as in (4) with triggering conditions given by 

Ti(xi,ei,di) =mm{T il (xi,ei),Ti2(x i ,ei,di)} , (45) 

for all i € Af. Then the origin is a globally uniformly asymptotically stable 
equilibrium for (44), if there are constants Kj > 0, j € M such that at the 
triggering times tk, which are implicitly defined by (44) and (45) as described 
in Section 3, the following condition is satisfied: 

\\\itfk)\\-dM)\<«A*M)\\ («) 

where dj (t k , ) is defined by (43) . In particular, no Zeno behavior occurs. 

Proof: Before we can use Theorem 3 respectively Theorem 6 to conclude 
stability, we have to exclude the occurrence of Zeno behavior. First note that 
condition (45) triggers an event if and only if Tn > and Tj2 > respectively 
condition (8) and (36) are violated. Now assume that the jth subsystem induces 
Zeno behavior. For simplicity, we omit the index j of the triggering times t 3 k . 
Hence, let tk the triggering times of the jth subsystem and t* = lim^oo tk the 
finite accumulation point. From Lemma 3 we know that the jth subsystem has 
to approach the equilibrium, i.e. lim t( ._> t » Xj(tk) = 0. Lemma 4 tells us that for 
all M there exists a k* such that for some k > k* 

> M\\ Xj (t k )\\ . (47) 

tk - tfc-i 

As discussed in the proof of Lemma 2, the full state x £ A C A. But the 
knowledge of x is not available to a single subsystem. Hence, we take £(tfc) € 
A(j, Xj(tk), dj(tk)) as in the definition of W as an approximation for the states 
of the other subsystems. For this £ we can deduce together with the Lipschitz 
continuity of 0j and (47) 



LjmpiUiWUxjitkm > QMil- ■ ■ - INI. ■ • • - II6HI) > 

^f- x / tk -^ kjWxj^W > (M - K3 )\\x 3 (t k )\\ ■ (48) 
And hence for the k given in (47) 

max{||&(t fc )H, \\xj(tk)\\} > ^f-^||^-(tfc)|| • (49) 

Now choose 

M > max{Kj + Lj,Kj + LjKj} , (50) 

where Kj is the Lipschitz constant of ipj 1 o fjj. From Lemma 4 we know that 
this choice of M yields a k* such that we can conclude together with (49) 
max^- { || ||, ||a;j(t fe )||} = max^j ||&(£fc)|| for some k > k* . For this k we 
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want to show that the corresponding tk is not a triggering time. 
To this end we use (48) and (47) to get 

max||&(i fe )|| > — (1 - — ) . (51) 

«7=J J-ij M tk — tfc-1 

Note that for the jth subsystem (51) is true for all £ € A and therefore by the 
definition of W 



Lj M tk — ifc-i 

Using the latter inequality and the Lipschitz constant for tpj 1 ofjj we can bound 
Tj2 by 

r j2 < Aj||ej(i fe )|| - — (1 - — ) . 

M tk — tk-i 

From the definition of ej(tk) — Xj(tk-i) — %j{tk) w e arrive at 

t < k (t \ „ ft mi 1 n K ^ H x jfa) -^fa-i)|| 
r j2 < Kj\\Xj(t k ) - a;j(ife_i)|| - —(1 - — ) _ - . 

We may assume that fc* is sufficiently large so that tk — tk-i < M~ x for all 
k > k*. Together with (50) we obtain 

Kj < Ljitk-tk-i)^ ~ 

and hence Tj 2 < in contradiction to the assumption that tk is a triggering 
time. Because the only further assumption on the solution of (44) and (45) is 
the occurrence of Zeno behavior, the aforementioned contradiction shows that 
Zeno behavior cannot occur. 
To conclude stability define 

I(x,e):={j€Af : Vj( Xj ) > Xj(\\ej\\)} , 

J(x,e):={jeAf : V(x) > %(||e,||)} , 

and 

J(x,e) := {j € N : ^(W&xj,^)) > %(|| e ,||)} . 

Note that the triggering condition Tj ensures that j G I U J . For j e / wc can 
use exactly the same reasoning as in Theorem 3. 

Lemma 2 tells us that from j £ I and j e J we can deduce j e J. For the 
case j G J wc can adopt nearly the same reasoning as in Theorem 3. Only 
the reasoning for the existence of a Lyapunov function for the overall system 
changes. In Theorem 3 it can be deduced from Theorem 2 whereas here we have 
to use Theorem 6 to conclude the existence of a Lyapunov function. The rest 
of the proof can be copied word by word from Theorem 3. This ends the proof. 

□ 

Remark 8 The advantage of parsimonious triggering is twofold. First it allows 
us to exclude the occurrence of Zeno behavior and second it may lead to fewer 
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transmissions compared to the triggering condition given in Theorem 3. Com- 
pared to Theorem 5 where the same goal is achieved by the notion of practical 
stability, here we still achieve asymptotic stability, but we have to place more 
technical assumptions on the involved class K. estimates. 

Note that the set of indices j £ J\f for which condition (8) holds is a subset of 
those for which (36) holds (in other words I(x,e) C J(x,e)). But we cannot 
check condition (36) locally. 

Because of the conservatism we introduce by using T i2 instead of (36), trigger- 
ing condition Ti still makes sense. 

In a practical implementation Tn should be checked first, before T i2 is checked, 
because of the conservatism of T i2 and the possible cumbersome calculation of 
W. 

Remark 9 One possible drawback of the triggering condition given in Theo- 
rem 7 is that the condition on the approximation dj as in (46) might be too 
demanding. First note, that if tk — tk-\ is sufficiently small, (46) trivially holds 
true, because dj is the difference quotient. As Xi approaches zero, it could hap- 
pen that the difference tk — tk-i does not decline fast enough to ensure that (46) 
holds. A way to overcome this issue would be to adjust condition T i2 in such a 
way that it always tries to trigger an event as soon as it cannot be guaranteed 
that the approximation dj satisfies (46). 

9 Conclusion 

We presented event-triggered sampling schemes for controlling interconnected 
systems. Each system in the interconnection decides when to send new infor- 
mation across the network independently of the other systems. This decision 
is based only on each system's own state and a given Lyapunov function. Sta- 
bility of the interconnected system is inferred by the application of a nonlinear 
small-gain condition. The feasibility of our approach is presented with the help 
of numerical simulations. To prevent the accumulation of the sampling times in 
finite time, we propose two variations of the event-triggered sampling-scheme. 
The first is based on the notion of input-to-state practical stability, whereas the 
second compares the local error to an approximation of the Lyapunov function 
of the overall system to guarantee stability of the interconnected system. 

A Appendix 

A.l Proof of Lemma 2 

Proof: For later use define 

A{ 3 ,x 3 ,e 3 ,ij) = : Bee R ,H s.t. f j (^,g j {^ x '+^)) = ±j and 

V{^)>UM\)^^3] ■ (52) 

The set A(j,xj,ej,Xj) describes the set of all ^' Xj for which a pair (£ J ' X 3 ', e J ' ej ') 
exists that fulfills the right hand side of the jth subsystem for a given ij , Xj , ej 
and for which V(^' Xj ) > %(||ej||) for all i ^ j hold. As the system's state 
satisfies the dynamics, it holds that x e A. 
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Before we proceed, we want to show that A(j, Xj, ej, Xj) C A(j, Xj, dj). To this 
end take a £ e A(j,Xj,ej,Xj). Hence we have fj(^' Xi ,gj(^' Xi + e J ' ej )) = x j- 
Taking the norm and using (37) yields 

ej(Ui\\,...,Un\\) > \\m j ' Xi ,9j(^+e j ' x m\ = > dj-iijWxjW , 

where the last inequality follows from the condition on the approximation for 
|| Xj || . And we can conclude A (j, xj , ej ,Xj) C A(j, Xj ,dj). 
From condition (39) we can deduce 

^J 1 °»KIMI) < W(j,Xj,dj) = min{max||&|| : £ € A(j,Xj,dj)} < 

minjmax H^ll : £ £ A(j,Xj,ej,Xj)} < max.\\xi\\ . (53) 

The second inequality follows from A(j,Xj,ej,Xj) C A(j,Xj,dj) and the last 
can be deduced from x € A. Now we can rewrite (53) to get 

%(IN|) < ^(maxll^H) . 

^3 

With the help of (39) and the definition of tpj we arrive at 

^(Vjfa)) < %(ll e jll) < max^r 1 ° aii(lkill) < maxcr^ 1 ^^)) , 

1^3 

where the last inequality follows from Assumption 1. Considering the first and 
the last term in the chain of inequalities above it is easy to see that the jth 
subsystem does not contribute to the Lyapunov function of the overall system 
and we conclude max,^ , a^ 1 (Vi(x i )) = ma,x ie ^f a^ 1 (Vi(xi)) = V{x) and the 
proof is complete. □ 



A. 2 Proof of Lemma 3 

Proof: Denote t* = lim/^oo tk ■ By definition of the triggering condition 
we have for each k an index i{k) € Af such that 

Vi(k){xi( k ){tk)) = Xi(fe)(l|ei(fe)(ife)||) • 

Choose i* E Af such that i{k) = i* for infinitely many k. Such a i* exists 
because Af is finite and k ranges over all of N. Let K be the set of indices for 
which i(k) — i* . For ease of notation let K = {s\,S2, . . . }. By Theorem 2 V is 
a Lyapunov function for the event triggered system on the interval [0, t*). Thus 
the trajectory £|[o,t») is bounded and e| [o,t*) is bounded because Xi(ll e i(^)ll) < 
Vi(xi(t)) for alii e Af, t € [0,t*). It follows that Uj|[o,t*)) is bounded and so ±i 
is bounded on [0,t*) for all i € Af. 

Then we have by uniform continuity of X{* on [0, t*) that the following limit 
exists 

lim xA\Ms k )\\) = lim Vi.(xi.(« fc )) = MM**)). ( 54 ) 

k — Voo k— >oo 

By definition e f . (s£) = 0. By (4) we have that = — ij* almost everywhere on 
(sfe, Sfe+i). Since ±j* is bounded and s^+i — Sk — > 0, then condition ej* (s^ ) = 
implies that 

ei*(T)dr = / ei.(r)dr 

which goes to for fc — > oo. Hence by (54) we obtain that Vj* (xi» (t*)) = 0. 
This shows the assertion. □ 
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A. 3 Proof of Lemma 4 

Proof: The proof will be by contradiction. To this end assume that for 
some fixed M > and all k sufficiently large we have 

\\xi{t k +i) - Xi(t k )\\ < M(t k +i - t k )\\xi(t k +i)\\ ■ (55) 

The evolution of Xi between ti and t k for k > I can be bounded by using a 
telescoping sum, the triangle inequality applying (55), and a judicious addition 
of 0: 



\\Xi(t k ) - Xiit^W < 
fc-i 

3=1+1 

If we choose D > and a k' such that < , — , 4lt 1 -. T < D for all k > k' we 

can rewrite the latter to 

fe-i 



\\xi{t k ) - Xi{tt)\\ < 1 _ M{ l_ tk i) E M(t j -t j - 1 )\\x i (t j )-x i (t l )\\ + 

xmw . 



3=1+1 

M(t k -tt) 



1 - M(t fc - t fc _i) 

Using the discrete Gronwall inequality (see e.g., Theorem 4.1.1 from [1]) yields 

wxi(t k ) - xmw < 1 

1 - iW(t fc — tfc-ij 

Exploiting that f& is a monotone sequence, that < 1 _ M ^ t 1 k _ tk ^ < and that 
1 + x < e x for all x e R and collapsing the telescoping sum again gives 

WxiW-XiWW < (MD(t k - t,) + M 2 D 2 (t k -i - tl fe MD ^-^) |M* ; )|| • 

s v ' 

-c 

Because of the finite accumulation point t* , there exists an k* > k' such that 

\\xi(t k ) - Xi (ti)\\ <C\\ Xi (ti)\\ 

for all k > I > k* with C < 1. Realizing that this contradicts ||a;j(i*)|| = 
finishes the proof. □ 
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